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Abstract 

We study SU{3) flavor-breaking corrections to the relation between the octet baryon masses and 
the nucleon-meson CP-violating interactions induced by the QCD 6 term. We work within the 
framework of SU{3) chiral perturbation theory and work through next-to-next-to-leading order in 
the SU{3) chiral expansion, which is At lowest order, the CP-odd couplings induced by the 

QCD 0 term are determined by mass splittings of the baryon octet, the classic result of Crewther et 
al. We show that for each isospin-invariant CP-violating nucleon-meson interaction there exists one 
relation which is respected by loop corrections up to the order we work, while other leading-order 
relations are violated. With these relations we extract a precise value of the pion-nucleon coupling 
go by using recent lattice QCD evaluations of the proton-neutron mass splitting. In addition, we 
derive semi-precise values for CP-violating coupling constants between heavier mesons and nucleons 
with ~ 30% uncertainty and discuss their phenomenological impact on electric dipole moments of 
nucleons and nuclei. 


1 



I. INTRODUCTION 


Violation of time reversal (T), or, equivalently, violation of the product of charge conju¬ 
gation and parity {CP) is one of the ingredients pQ needed to explain the matter-antimatter 
asymmetry of the visible universe. The Standard Model (SM) of particle physics contains two 
sources of CP violation, the phase of the CKM matrix and the QCD 9 term. The phase of 
the CKM matrix explains the observed CP violation in K and B decays |2], but appears to 
be too small for electroweak baryogenesis (see j3] and references therein). The QCD 9 term 
is severely constrained by the non-observation of the neutron electric dipole moment (EDM). 
The current limit on the neutron EDM, \dn\ < 2.9 ■ 10“^^ e fm |1], bounds 9 to be small, 
9 < 10“^°, the well known strong CP problem. 

The viability of attractive, low-scale baryogenesis mechanisms such as electroweak baryo¬ 
genesis thus requires new sources of CP violation. With the assumption that new physics is 
heavier than the electroweak scale and that there are no new light degrees of freedom, new 
sources of CP violation appear as higher-dimensional operators in the SM Lagrangian, sup¬ 
pressed by powers of the scale at which T violation arises. These operators involve SM 
particles and respect the SM gauge symmetry. In the quark sector, they are at least dimension 
six, and are suppressed by two powers of [3 E]. EDMs of the nucleon, light nuclei, atomic, 
and molecular systems are extremely sensitive probes of such flavor-diagonal T-violating {P^ 
operators, for reviews see Refs. [7H9]. The current generation of experiments probes scales of 
TeV (or higher), and provides powerful constraints on new physics models, complementary to 
direct searches of new physics at colliders. Furthermore, a vigorous experimental program is 
under way m to improve current bounds on the neutron EDM [1] by two orders of magni¬ 
tudes, to measure EDMs of light nuclei at the same or even higher level of accuracy mm, 
and to improve the bounds on EDMs of diamagnetic atoms, like ^^®Hg [I3], ^^®Xe [ll] and 

225Ra m- 

The extraction of robust information on possible new sources of CP violation from EDM 
measurements involves dynamics on a large variety of scales, from the new physics scale Mj^ 
to the electroweak (EW) and QCD scales, down to the atomic scale. The step to hadronic 
and nuclear scales involves nonperturbative strong matrix elements which are often poorly 
known leading to large uncertainties [H]. In recent years, lattice QCD has made progress in 
computing the nucleon EDM induced by the QCD 9 term [TBHTB] . while the study of the 
nucleon EDM generated by higher-dimensional operators is still in its infancy (an exception is 
the quark EDM [I9j). Another important ingredient for the study of EDMs are^ pion-nucleon 
couplings, which determine the leading non-analytic contribution to the nucleon EDM m- 
In addition, they generate ^ long-range nucleon-nucleon potentials, contributing to EDMs of 
light nuclei |21l [22] , and diamagnetic atoms [Hi- 

The direct calculation of ^ pion-nucleon couplings on the lattice, both from 6 and 
dimension-six operators, is difficult. Some information can be gained by the study of the 
momentum dependence of the electric dipole form factor (EDFF) [231121], but the most recent 
lattice calculations are performed at too large momenta for a reliable extraction [ISICSII2SI. 
Fortunately, in some cases other methods exist to extract the values of the pion-nucleon in¬ 
teractions. For chiral-symmetry-breaking sources, like the QCD 9 term, or the light-quark 
chromo-electric dipole moments (qCEDMs), the pion-nucleon couplings are intimately related 
to UP-even, chiral-symmetry-breaking effects. In the case of the QCD 9 term this was realized 
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in Ref. ra. that expressed the isoscalar 'J' coupling go in terms of mass splittings of the octet 
baryons. In particular, in SU{2) Chiral Perturbation Theory (yPT) it is possible to relate go 
to the neutron-proton mass difference induced by the quark mass difference — m„, which 
we denote by 6mi^. All the information on non-perturbative dynamics entering go can thus be 
extracted by computing a property of the baryon spectrum, the nucleon mass splitting, a task 
for which lattice QCD is particularly well suited. Indeed, existing calculations of the nucleon 
mass splitting allow a determination of go with 10% accuracy [26l EZ], if one considers only 
lattice uncertainties. Similar relations between ^ pion-nucleon couplings and modihcations to 
the meson and baryon spectrum can be derived in the case of the qCEDM operators [3 EH], 
and provide a viable route to improve the determination of couplings that, at the moment, 
are only known at the order-of-magnitude level [TJ [H EHj. 

The relations between ^ couplings and baryon masses strictly hold at leading order (LO) 
in yPT. Furthermore, if one considers the strange quark as light, and extends the chiral group 
to SU{3) X SU{3), more LO relations can be written, e.g. go can be expressed in terms of 
the mass difference of the S and S baryons. Using the two LO relations leads to values of go 
that differ by about 50%, well beyond the lattice QCD uncertainty. Is this large difference due 
to an inherent uncertainty in the relation between the spectrum and the ^ couplings? Does 
this imply that the relations to the baryon spectrum can only be used for order-of-magnitude 
estimates of the T couplings? 

In this paper we investigate these questions and seek to quantify the SU{3) flavor-breaking 
corrections between the baryon masses and^ couplings induced by the QCD 6 term. We work 
in the framework of SU{3) heavy baryon yPT [HUl ET] and compute higher-order corrections 
in the chiral expansion. We show that most LO relations are badly violated, already at 
next-to-leading order (NLO) and cannot be used for reliable extractions of the meson-nucleon 
couplings. However, for all isospin-invariant ^ couplings there exists exactly one relation that 
is preserved by all loop corrections up to next-to-next leading order (N^LO). By using the 
relations that are not violated by SU (3) flavor breaking, a precise extraction of the couplings 
is possible irrespective of the convergence of SU{3) xPT through this order. In the case of 
go, the preserved relation is to ^m^r, while the relation to the mass difference of the S and 
S baryons receives large NLO and N^LO corrections, which show little sign of convergence. 
Expressing go in terms of the S and S masses overestimates the coupling by about 50%, well 
outside the uncertainty which is determined with 

For isospin-breaking couplings, such as the isovector pion-nucleon coupling gi, we were not 
able to identify any robust relation that does not receive large violations already at NLO. We 
are forced to conclude that SU{3) heavy baryon xPT does not provide a reliable method to 
extract this important coupling from known matrix elements. 

In this work we only focus on CP violation from the SM QCD 6 term, leaving higher¬ 
dimensional operators arising from possible BSM physics for future work. However, our results 
are also relevant for scenarios of BSM physics where the strong CP problem is solved by a 
Peccei-Quinn mechanism. In this case, an effective 6 term can be induced proportional to 
any appearing higher-dimensional CP-odd sources [7l[32|. Other BSM scenarios involve cases 
where parity is assumed to be an exact symmetry at high energies, requiring 6 = 0 [33], 
while a calculable contribution to 6 is induced at lower energies once parity is spontaneously 
broken, see for instance Ref. [3l]. In any case, a quantitive understanding of the low-energy 
consequences of the 6 term is necessary to unravel the underlying source of CP violation once 
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a nonzero EDM is measnred [35] and to test scenarios, snch as the one in Ref. [36|, where a 
small bnt nonzero 6 term is expected. Our values of the OP-odd pion-nucleon couplings can 
also be used for more precise limits on axion searches [371 [38] . 

This paper is organized as follows. In Section we review the closely related chiral sym¬ 
metry breaking and ^ sectors of the xPT Lagrangian. In Section III we discuss baryon masses 
and^ couplings at tree level, and identify the relations between masses and couplings imposed 
by SU{3) symmetry. In Section IV we study NLO corrections to masses and^ couplings, and 
identify which relations are respected by NLO loop corrections. In Section jV Aj we discuss in 
detail N^LO corrections to the nucleon mass splitting, including, for the first time, decuplet 
corrections. The expressions of N^LO corrections to the mass splittings of the S and S baryons 
and to the octet baryon average masses are relegated to Appendices and [^ In Section jVBj 
we examine N^LO corrections to Qq, and show that all the loops at this order are related to 
contributions to the nucleon mass splitting. In Sections V C and V D we discuss the remaining 
^ nucleon couplings. In Section VI we use the conserved relations to determined the value of 
the couplings induced by the QCD 9 term, and discuss the impact of our analysis on the 
nucleon electric dipole form factor and on the ^ nucleon-nucleon potential. We conclude in 
Section IVIIl 


II. QCD AND EFT LAGRANGIAN 

At the QCD scale, p ~ 1 GeV, heavy gauge bosons, the Higgs and the heavy quarks can be 
integrated out, and the SM Lagrangian involves gluons, photons, and three flavors of quarks. 

£qcd = . (1) 

g is a triplet of quark helds q = {u, d, s), and are the photon and gluon held strengths, 
and is the SU{3)c x U{l)em covariant derivative. The hrst three terms in Eq. ([^ are the 
photon, gluon, and quark kinetic terms. Without loss of generality the quark mass matrix can 
be expressed in terms of a real, diagonal matrix Ai = diag(m„, m^, mg), and a common phase 
p. The last term in Eq. Q is the QCD 9 term. Despite being a total derivative, it contributes 
to physical observables through extended held conhgurations, the instantons [39]. The two 
CF-violating parameters in Eq. ([^, 9 and the phase p, are not independent, and CP violation 
is proportional to the combination 9 = 9 — Ufp, where rij = 3 is the number of havors of light 
quarks. This can be explicitly seen by performing an anomalous Ua{^) axial rotation. With 
an appropriate choice of phase, the 9 term can be completely eliminated, in favor of a complex 
mass term. The residual freedom of performing non-anomalous SU{3)a axial rotations can be 
used to align the vacuum in presence of CP violation to the original vacuum of the theory. If 
the complex mass term is the only^ operator in the theory, vacuum alignment is accomplished 
by making the complex mass term isoscalar m. At the level of the meson Lagrangian, the 
condition of vacuum alignment is equivalent to setting the leading-order coupling of the pion 
and r] meson to the vacuum to zero. 

After vacuum alignment, the QCD Lagrangian in the presence of the 9 term reads 

Gqcd = -jiv.-F'"' - jG^^G”"" + iiipii q , (2) 
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where we denote 


mumdUis __ m(l - e^) 

ms{mu + md)+mumd 2+^(l-£2) 

with 2m = rriu + md and £ = (m^ — mu)/{md + mu)- When providing numerical results, we 
take the values of these quantities from the most recent lattice average by FLAG (quoted in 
MS scheme at /i = 2 GeV when relevant) [H]: 

Tfl 

m = 3.42 ± 0.09 MeV , — = 27.46 ± 0.44 , e = 0.37 ±0.03 . (4) 

m 

The QGD Lagrangian is approximately invariant under the global chiral group SU{3)l x 
SU{3)r. Ghiral symmetry and its spontaneous breaking to the vector subgroup SU{3)v lead 
to the emergence of an octet of pseudo Nambu-Goldstone (pNG) bosons, the pion, kaon, 
and Tj mesons, whose interactions are dictated by chiral symmetry. The quark mass and the 
QGD 6 term break chiral symmetry explicitly. Ghiral invariance can be formally recovered by 
assigning the mass term the transformation properties 

A4 + im^9 —)■ R(Ai + , M. — im^,9 —>■ L{M. — im^9)R^ , (5) 

under a SU{3)l x SU{3)ji rotation. The QGD 9 term thus induces;?^ interactions between 
pNG bosons, and pNG bosons and matter fields, that can be constructed using the same 
spurion fields employed in the construction of the meson and baryon mass terms. We refer 
to Refs. |2ni [2H1112H1S] for more details. In the next Sections we give the meson and baryon 
yPT Lagrangians relevant to the calculation of;2^ baryon-pNG couplings at N^LO. 


A. Meson Sector 

The constraints imposed by chiral symmetry and its spontaneous and explicit breaking on 
the interactions of pNG bosons can be formulated in an effective Lagrangian, yPT |121 H5L 
[49] . whose construction is well known. We adopt here the notation of Ref. [45]. In the 
absence of explicit chiral symmetry breaking, the interactions of pNG bosons are proportional 
to their momentum, q, which guarantees that low-momentum observables can be computed 
as a perturbative expansion in where A^ is a typical hadronic scale, A^ ~ 1 GeV. 

The quark masses explicitly break chiral symmetry, giving masses to the pNG bosons and 
inducing non-derivative couplings. However, the breaking is small and can be incorporated in 
the expansion by counting each insertion of the quark mass as q‘^. 

We assign each term in the yPT Lagrangian an integer index, that counts the powers of 
momentum or of the quark mass. The LO meson Lagrangian contains two derivatives or one 
light quark mass insertion, and is given by 

Ci^) = ^TT{d,Ud>^U^) + ^Ti[U^X + Ux^] , ( 6 ) 

where Fq is the pion decay constant in the chiral limit. Beyond LO, SU{3) breaking corrections 
break the degeneracy of the pion, kaon, and rj decay constants. We denote by Ft, and Fk the 
empirical pion and kaon decay constant, ±4 = 92.2 MeV and Fr = 113 MeV [2|. can be 
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expressed in terms of Fk and and we use = 
unitary matrix 

f/(7r) = u{7iY = exp 


l.SFjr |12]. In Eq. 
2i7r\ 

’ 


([^ we introduced the 
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where tt are the pNG boson fields 


TT = TT t,. 


1 
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00 




( 8 ) 


(t°‘)ij are the generators of SU{3), a = is the octet index and i,j = are 

indices of the fundamental representation of 517(3). Under a SU{3)l x SU{3)r transformation, 
the pNG held has a complicated non-linear transformation, while U transforms simply as 
U —)■ RUL"^. The hrst term in Eq. (§ is chirally invariant. The second term, with y = 
2B{M. + im^6), is the realization of the quark mass term which, with the transformation 
properties in Eq. ([^, is also formally invariant. 

Eq. ([^ induces the leading contribution to the pion, kaon, and rj meson masses. 

m‘t± = 2Bm , 

TT-^ * 

^ 2^2 

mlo = 2Bm - B -^ , 

rus — m 

n 2 , _, ffi^e^ 

m = -B{2ms + m) + B -^ , 

' 3 rris — m 

m']^± = Birris + fh — me) , 

m^o = B{ms + fh + me) . (9) 


When working in the isospin limit, we will denote = B{ms + rh). At the order we are 
working, we need the meson masses only at LO, and, for numerical evaluations, we will use 
the PDG values mjr± = 139.6 MeV, mK+ = 493.7 MeV, mxo = 497.6 MeV, = 547.9 
MeV [2] . The experimental error on the meson masses is always negligible compared to other 
uncertainties in the calculations, and we can neglect it. 

The relation between the physical tto and rj and the pNG bosons tts and vrg is determined, 
at LO, by the n - r] mixing angle 0, 


TTo = COS 0 TTs -7 siu 0 TTg , 

7] = — sin 0 TTg -f COS 0 TTg , (10) 


with 


me 


( 11 ) 


^/3 2{ms - m) 

Beyond lowest order, rj - n mixing cannot simply be described by a mixing angle 

In the sector, vacuum alignment eliminates pion and rj tadpoles in LO. In SU{2) yPT, 
vacuum alignment eliminates all LO three-pion vertices. On the other hand, in SU{3) yPT 
the meson mass term induces a three-pNG vertex of the form 


r = 


B 


m^9 (f^^TTa'n'bT^c , 


3Fr 


( 12 ) 


6 








where are the constants determined by the anticommntator of SU{3) generators 
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(13) 


The interaction in Eq. (12) involves one r] and two pions, one i] and two kaons, or one pion 
and two kaons, and indnces the CP-odd decay rj —)■ tttt. Limits on this branching ratio allow 
to pnt a bonnd on 6, thongh several orders of magnitndes less stringent than the bonnd from 
the nentron EDM [5T]. Three-pion interactions also arise at LO, bnt they are proportional to 
the rj - 7T mixing angle, and vanish for large rUg. 

The meson Lagrangian is well known |12]. At the order we are working, we only 

need the terms 


4^) D L 4 Tr(a^pt Tr (x^U + U^x) + ^sTr {d^U^ {x^U + U^x)) 

+Lj (Tr(Px^ -XU^)Y + LsTt {Ux^Ux^ + xU^U^) • (14) 

L 4 and L 5 contribnte to the pNG wave fnnction renormalization and to the renormalization 
of pQ. L^ and Lg generate pion and rj tadpoles, that contribnte to T pion-nncleon conplings 
at N^LO. For ^ baryon-pNG conplings, the dependence on L 4 and L 5 cancels between the 
wave fnnction renormalization and the corrections to Pq. Lj and Lg have been determined 
from global fits to meson data [521 E3]. We nse the NLO fits in Ref. [53], which give L 7 = 
(—0.3 ± 0.2) ■ 10“^ and Lg = (0.5 ± 0.2) ■ 10“^. Lg is scale dependent, and it is evalnated at 
the scale /i = 770 MeV. 


B. Baryon Sector 

The inclnsion of baryons in yPT has been derived in a large nnmber of papers, for instance 
Refs. [301 EH l5H - f57j . The baryon octet can be inclnded in a way consistent with the chiral 
expansion by working in the non-relativistic limit and removing the large, inert octet mass 
tiib [301[3T]- The mass splittings of octet states vanish in the chiral limit, and scale as 
yPT can be extended to inclnde the decnplet baryons at the price of introdncing a new scale 
A, the decnplet-octet splitting, which does not vanish in the chiral limit nor can it be rotated 
away [sni EZj. This octet-decnplet splitting scales as l/A), in the large Nc expansion [35L 
EO]. The explicit inclusion of the decnplet is necessary for the chiral expansion to respect 
the 1/Nc counting rules [ 6 l], and the predictions from a combined SU{3)-l/Nc expansion are 
phenomenologically well satisfied in lattice QGD calculations [62] ES]- 

In the heavy baryon formalism, the lowest-order chiral-invariant octet and decnplet baryon 
Lagrangian is given by 

£(1) = Tr {iBv ■ VB) + P Tr {BS^ 4, P]) + P Tr {BS^ {u^, B}) 

-iTv ■ VT^ + ATt^ + ^ {Tu^B + Bu^T^) + UTS^u.T^ , (15) 

where Vfj, and denote the heavy baryon velocity and spin. D and P are the octet axial 
couplings, D ~ 0.8 and P ~ 0.45, and D + F = qa = 1.27. C is the decnplet-octet axial 
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coupling, C ~ 1.5. H is the decup let axial coupling, which does not play a role in our 
discussion. The matrix B denotes the octet baryon held 


B 




1 vO 




._Lyo _i_ _LA 



( 16 ) 


The decup let held carries three completely symmetrized fundamental indices. The pNG 

bosons appear through the combinations = uHdfj,u — uid^u^ and, in the covariant deriva¬ 
tives, = \{ud^u^ + u^d^u). The chiral covariant derivatives are 

V^B = d^B + [V^,B] , (17) 

wr 6 , = + (K),, Tt;, + (K,)„ TA + (V.)u % . ( 18 ) 

Invariant terms involving the octet and decuplet baryons are constructed using the contractions 


, T^u.B = . ( 19 ) 

Baryon mass terms and ^ couplings appear in the SU{3) Lagrangian at 0{q‘^), and they 
are given by [SU] 

= boTT{BB)TTx+ + bnTT{B{x+,B})+bFTT{B[x+,B]) 

+bc Tx+T, + 6a Tr(xl7t + x^U) T% , (20) 

where x+ = u^X'^^ + ux^u. bo and 6 a denote common shifts to all octet and decuplet masses, 
and do not give rise to nucleon couplings with one pNG boson, bo, bp, and be induce 
splittings between the diherent octet and decuplet states, and give rise to ^ baryon-pNG 
interactions. The LEGs bo, bp, and bp scale as A“^, and do not depend on the quark masses. 
When including decuplet corrections bo, bp, and bp must be interpreted as series expansions 
in the octet-decuplet splitting A [MHOTj 




1 




A 



X 



( 21 ) 


where the higher orders in A/A^^ arise from hnite contributions and are needed to absorb 
the divergences arising from diagrams with decuplet intermediate states. There is no sense 
in keeping track of the hnite A dependence in the LEGs as these corrections are quark mass 
independent, and thus not discernible with present lattice QGD calculations at hxed N^. = 3. 

Baryon mass splittings and^ couplings receive 0{q^) corrections from one-loop diagrams 
involving the chiral-invariant interactions in Eq. (15), and chiral-breaking interactions from 
Eq. (|^. This is distinct from SU{2) xPT in which the 0{q^) corrections cancel in the isospin 
mass splitting. At 0{q'^), one has to consider one-loop diagrams involving operators in the 
0(q‘^) Lagrangian, and tree-level diagrams with one insertion of the 0{q‘^) Lagrangian. Besides 
the mass terms in Eq. (20), contains relativistic corrections to the interactions in Eq. 
(15), and baryon-pNG interactions with two derivatives. The operators are listed in Ref. [68] 


8 





and here we give only those relevant to our discussion. There are four operators containing 
two derivatives of pNG helds. 


= 6 iTr [B K, K, i?]]) + &2Tr {B 5}]) + fegTr {B 5}}) 

+hf{Ii[BB]Ti[u^u^^] . ( 22 ) 

The effects on the baryon masses of operators similar to 6i, ..., feg, but with replaced by 
V ■ u, can be accounted for by a redehnition of and of the 0{q^) LECs. For this reason, we 
do not include these operators explicitly. The relativistic corrections are 

£2 = -^Tr (BS„ [iD^ {v . n, B}]) - ^Tr (BS„ \v ■ u, B]]) 

-^IV (BS, [v ■ u, liB", B|]) - ^IV (BS, {v ■ u, [iD^, B]}) 

Jj2 _ 3^2 _ Jj2 _ 

H——-Tr (B [n • u, [n • u, B]]) — — - Tt{BB)Tt{v ■ uv ■ u) 

24:171 Q lIttiq 

I®"- I*"' 

DF 

-^^Tr {n-M,^}]) . (23) 

We hnd that the contribution of recoil corrections to D and F to the baryon masses are small. 
The relativistic corrections to the octet-decuplet coupling C can be removed using the LO 
equations of motion |U3] - IU7] . 

The 0{q^) Lagrangian is 

= diTr {B [x+, [x+, 5]]) + {B [x+, {x+, B}]) + dgTr {B {x+, {x+, B}}) 
+d^TT{Bx+)Tr{x+B) + dgTr {B [x+, B]) Tr(x+) + deTr {B {x+, B}) Tr(x+) 
+d7Tr(Si?)Tr(x+)Tr(x+) + dsTT{BB) Tr(x^) 

TdgTr {B [x_, [x-,B]]) + dioTr {B [x_, {x-, 5}]) 

+diiTr {B {x-, {x-, B}}) + diaTr (i?X-) Tr (x-B) 

TdisTr {B [x-, B]) Tr(x-) + duTr {B {x-, B}) Tr(x-) 

+di5Tr(SS)Tr(x-)Tr(x-) + di6Tr(SS)Tr(x!) • (24) 

di,... ,ds were constructed in Ref. [6H], and contribute to baryon masses and splittings. The 
operators dg,..., die involve two insertions of X- = 7xu^ ~ ux7. The GP-even parts of 
these operators do not contribute to baryon masses and mass splittings, but do contribute to 
pion-nucleon scattering. The CP-odd components give 0{q‘^) corrections to 7 baryon-pNG 
couplings. 


III. OCTET BARYON MASSES AND 7 COUPLINGS AT TREE LEVEL 


The Lagrangian (20) realizes the leading effects of the light quark masses in the baryon 
sector. The light quark masses induce splittings between the octet and decuplet states, and. 
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in the presence of the QCD 6 term, cause the appearance of^ couplings between baryon and 
pNG bosons. The LO corrections to the baryon masses are well known (see, for example. Refs. 
[6^ 170]). we give them here in order to make the connection with^ couplings explicit. 

The nucleon, S and S mass splittings are given by 

= rrin — mp =—SBrhe {bp + Bd) , 

= m=- — m=o = —8Bfhe {bp — bp) , 

= ms+ — rnp,- = IGBrhebp , (25) 


where we introduced the superscript (0) to denote that these are the leading contributions. The 
three mass splittings are not independent, but are related by the Coleman-Glashow relation 

[ZH 

= 0 . (26) 

Neglecting 0{e‘^) corrections to the and A masses, the isospin-averaged masses of the 
nucleon, S, A, and S baryon are 


A^^^rriN 

A^°Wa 


rrin + rup 


- mp = -4:B{ms{bo + bp - bp) + m{2bo + bp + bp)) , 


rri’z- -|- m^o 


mp = -4:B{ms{bo + bp + bp) m(26o + bp- bp)) , 


ms+ + m-po + ms- 


mp =—4:B{{ms + 2m)bo + 2mbp) , 


mA - mp = -AB ( (m* -h 2m)bo + - {m + 2ms)bp 


(27) 


where all the masses are measured with respect to mp, the common octet mass in the chiral 
limit. Finally, mp gets a corrections proportional to ms + 2m 


^ 2A*^°)mAr -f- 2A^^'>m'= -|- 3A*^°Ws + A^^Wa 

A} 'mp =-=- 


—AB{ms + 2m) 




(28) 


In the presence of a 0 term, the operators in Eq. (20) induce^baryon-pNG couplings. The 
couplings of the greatest phenomenological interest are pion-nucleon couplings. Besides giving 
a LO contribution to the nucleon EDM, the isoscalar non-derivative pion-nucleon coupling qq 
induces a.T nucleon-nucleon potential, which is expected to give a sizeable, when not dominant, 
contribution to EDMs of light nuclei with N ^ Z [211 [72l [73] . Furthermore, EDMs of heavier 
systems, like are commonly computed in terms of three non-derivative pion-nucleon 

couplings [H|. 

Introducing the nucleon doublet, N = (pn)^, we write the^ pion-nucleon couplings as 


C. = - ^Nr ■ 77A - ^ttoNN - + ... , (29) 

where r are the Pauli matrices, and ... include terms with more derivatives. Notice that we 
defined the couplings in terms of the physical pion decay constant F^, rather than Fq. The 
difference between iA- and Fq is an N^LO correction. 
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At tree level, the T pion-nucleon couplings are expressed in terms of the LECs bo and 
bp- bo does not generate tree-level T couplings with only one pNG boson, but does induce 
couplings with at least three pNG, which are relevant at one loop, go, gi, and g 2 are given by 


-(0) 

% 

= -3B{bD + bp)m^,6 , 

(30) 

-(0) 

91 

= 3B{bpi-3bp)-^mJ) , 

V 3 

(31) 

-(0) 

92 

= AB{bp) + bp)(tP‘m,^9 , 

(32) 


where 0 is the g - tt mixing angle dehned in Eq. (11). From Eqs. (25) and (30), we see that 


go is related to the tree-level contribution to the nucleon mass splitting. 




me 


2e 


m 


For mg 3> m, this is the same relation that holds in SU{2) 
one can also write [20] 


m.o 


We will show that both Eq. (33) and (34) are violated at N^LO. However, Eq. (33) is only 


m 

2l 


(33) 

In SU{3) yPT, at tree level 

(34) 


violated by hnite terms and by new LEGs appearing at while it is respected by all loop 

diagrams. On the other hand, Eq. (34) is already violated at NLO and receives much larger 
corrections. 

An isoscalar operator like 6 can generate the isospin-breaking couplings gi and ^2 only in 
the presence of some source of isospin violation. In SU{2) yPT this implies that gi and ^2 
are suppressed and appear at 0{q‘^) and 0{q^), respectively. In 517(3) xPT, the g - tt mixing 
angle (j) appears at LO which means that gi ^2 are formally LO as well. However, numerically 
they are suppressed by powers of fh/mg ~ 0.04. 

The coupling gi is particularly important for EDMs of nuclei with N = Z such as the 
deuteron. At LO, the combination of LEGs bp, — 3bp can be expressed in terms of baryon- 
mass splittings as 


-( 0 ) 

9i' = 


^/3me 
V3 [rrig - 


m] 


d 


d 


dm. 


2dm/ ' 


(35) 


where the masses and splittings are given in Eq. (27). Eqs. (31) and (35) were used in Ref. mi 
to estimate gi and its contribution to the deuteron EDM. Gontributions t o gi that are not 
suppressed by ffi/mg only appear at N^LO. However, as we discuss in Sect. V G, they can be 


as large as Eq. (35) 


In the case of the QGD 6 term, the coupling g 2 is suppressed by m^/mg, and because it is 
of little phenomenological consequence, we neglect ^2 henceforth. 
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Non-analytic LO contributions to the nucleon EDM in SU (3) yPT involve other ^ nucleon 
couplings [23]. Introducing the isospin doublet K = and an isospin triplet S = 

(S’*', E“), we can write the isospin-invariant^couplings between the nucleon and r] meson, 

and the nucleon, kaon and S or A baryon as 


C 


-^r]NN - ■ SA: - 

2Fri 2 Fk 


+ h.c. 


(36) 


Isospin-breaking couplings also arise at LO, but are suppressed by me/nis 
them in the following. At LO, the couplings in Eq. (36) are given by 


and we neglect 


-(0) 

= m^9 , 

Vs 

(37) 

-(0) 

9onek 

= 8B{bp - bD)mJ) , 

(38) 

-(0) 

9o NAK 

73 

(39) 


These ^ couplings can be expressed in terms of various combinations of baryon masses. We 
list some of them 



-( 0 ) 

9o NT,K 


V39<”' 


ONAK 


me 

-((i-4)) A'”’™") 






me 


mg — m 

me 

-(A(°)mA-A(°W)-^^ 

m. — m 


•> 




(40) 


(41) 


(42) 


where the approximate SU{3) symmetry enforces many other LO relations between corrections 
to the octet masses. The second equation in Eq. (40) is particularly interesting, since it relates 
Qor^ not directly to the nucleon mass, but to its derivatives with respect to and m. Through 
the Eeynman-Hellmann theorem, these derivatives can be related to the nucleon sigma terms 


O'Nq 


mq{N\qq\N) 


= mg 


5 Am AT 

dmq 


We can thus write 



I 

Y mg 


OnI \ 
2m j 


Am^9 , 


(43) 


(44) 
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FIG. 1: 0{q^) corrections to the baryon masses and baryon-pNG^ couplings. Plain, double, and 
dashed lines denote octet baryons, decuplet baryons, and pNG bosons, respectively. Dotted vertices 
denote CP-even couplings, i.e. the octet-pNG axial couplings D and F and the decuplet-octet-pNG 
coupling C. A square denotes a^ coupling. 


where am = ctnu + (^Nd- 

In the rest of the paper we show that in most cases the relations between baryon mass 
splittings and ^ baryon-pNG couplings break down already at NLO. The slow convergence 
of SU (3) baryon yPT then renders the usefulness of these relations to be qualitative only. 
However, for each isospin invariant coupling, there exist one relation that survives NLO and 
most of the N^LO corrections. These relations thus provide a powerful method to extract ^ 
couplings from well-known T-even matrix elements. The most important example is the link 
between go and Suin, which, as we discuss in Sections [IV] and |V} receives particularly small 
corrections. 


IV. OCTET BARYON MASSES AND 7" COUPLINGS AT NLO 

In Fig. [^we represent one-loop corrections to the baryon masses (diagrams [^a) and (b)), 
and to^ octet-pNG couplings (diagrams [^c) and (d)). In dimensional regularization, diagram 
l^a) is finite, and contributes to the octet masses and mass splittings at 0{q^). In particular, it 
affects the mass splittings Smisi, and Sm-^ through the kaon mass difference — = 

B{md — rriu), and the rj - tt mixing angle 0. At the same order, diagrams with intermediate 
decuplet states contribute via diagram [^b). These diagrams have UV poles that are linear in 
the octet-decuplet splitting A, and are absorbed by bp, and Bq. Diagram [^b) contributes 
to the nucleon mass splitting only through the kaon mass difference because of the vanishing 
of the octet-decuplet-T] axial coupling. The mass splittings of the S and S baryons receive 
contributions from both the kaon mass splitting and g - tt mixing. 
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All together we find 


^( 1 ) 


rriN = 




487rF2 


SttF^ 


73' 




1447r2F2 


(/(m^o, A) -/(mK+, A)) , 


(45) 




( A )2 + qdf - 3 F 2 ) 
487rF2 


(m' 


KO 


m 


3 . {D + ?,F){D-F) 0 ^3 3 . 


K+ j 


8nF^ 


:{m-m^) 




-liS^ A)) - 

(m%o - m^+) + — 


73' " 

(/(m^, A) - /(m^, A)) (46) 






247r2F2 


(/(mi^o, A) - /(mic+, A)) 


C 2 


(/(m^, A) - /(m^, A)). (47) 


127r2F2 43 

The loop function appearing in the decuplet diagrams is given by 

firriK, A) = A f-A^ + ^"^1-j + 4 fl2m^ - lOA^ + 3 {8m\ - 2A^) log 




m 


K 


+2{m\- — A^)^/^arccot- 


A 


7^1--^ 

L encodes the UV divergence, and is defined as 


(48) 


T = - + log 47 r - 7 e , 


(49) 


where 7 ^; is the Euler constant. For the spin projector for decuplet fields in d = 4 — 2^ dimen¬ 
sions, we used the definition of Ref. |H|. The poles are absorbed by defining the renormalized 
couplings Ifjj and 6 ^. We work in the MS scheme, and define 


bn = b 


'D 


A 




647r2F2 


L , 


= bF + A 


5C^ 


3847r2F2 


L . 


(50) 


In the limit A —>■ 0, the divergence disappears and / assumes the same form as the octet 
corrections 

hm^/(mx, A) = 7rm|- . (51) 


The NLO corrections in Eqs. (45), ( |46[ ), and ( |47[ ) that do not involve the decuplet agree with 
Ref. [7^. Both octet and decuplet corrections respect the Coleman-Glashow relation. 

In addition to the mass splittings, the tree-level relations between baryon masses and ^ 
couplings involve the octet isospin-averaged masses. Baryon masses in the isospin limit were 
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computed at NLO in Refs. [69l [70] and at N^LO in Refs. [65], IHH] IM] IZH] ES]. At NLO 


= ~ gQ^p2 (2(5A)^ - QDF + 9F^)m|- + 9{D + Fyml + {D - SF^m^) 




487r2F2 


(/(mx, A)+4/(m^, A)) , 


(52) 


A^^^ms = —TTT—^ (2(50“^ + QDF + 9F‘^)m\ + 9{D — F^ml + (D + 3F)^m^) 

yOTT ± n 


c' 


487r2F2 


(3/(mif, A) +/(m^, A) +/(m^, A)) , 


(53) 


A(^)ms = + F^)m% + A{D^+ QF^)ml + 4D^m^^) 






I4^4^^2p2 + 2/(m^, A) + 3/(m^, A))) 

((/}2 ^ ^ D^{3ml + mj)) 

(2/(mi^, A) + 3/(m^, A)) , 


(54) 


247rFo^ 




487r2F2 


(55) 


where the decuplet loop function / is given in Eq. (48). 

NLO corrections to gi are induced by the^three-pNG coupling in Eq. ( jl^ . This coupling 
is fixed at LO by the meson masses and does not involve a free coefficient. The relevant loop 
diagrams with octet and decuplet intermediate states are shown in Fig. Diagram [gc) is 
hnite, while [^d) is UV divergent. Both diagrams contribute to go, gi, and g 2 , although the 
last two couplings are suppressed by fh/ms and fri^/nil, respectively. The corrections to % 
and gi are given by 

-d) n n \ ~ 6DF - 3F^ m|.+ + mio + mK+rriKO 

g^ = Bm^O ■' 


247rF2 


niKO + niK+ 


{D - 3F){D + F) f ml + + ml^ ^ /(mA-o, A) -/(mx+, A) 


127rF2 


m^i + 


+ 


727r2F2 




(56) 


g[^^ = Bm^O 


5D^ - QDF + 9F^ 
327rF2 


mxo — mK+ + {m^o + mA'+) — 


g (D - SF^mr, - 5{D + F^m^ ^ {D‘^ + QDF - 3F^) m^ + mr,m^ + ml 


IQtiFq 


QttF^ 


mri + mj, 




327r2F2 


C2 


fimxo. A) - f'{mK+, A) + ^(/'(mA-o, A) + f\mK+, A))^ 
f{mr„A) - f{m^,A))' 


y3 67r2F2 


3/'(m^,A)-2^ 


mS — ml 


The function /' entering the decuplet corrections to gi is 


73 


(57) 


(58) 
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Notice that loops with only pions do not contribnte to qq at NLO, in accordance with the SU{2) 
resnlt of Ref. m- The piece proportional to {D + contribnting to is the same as 

fonnd in SU{2) yPT [22] once the LO identifications D + F = qa and (j)/y/S = {6ml)/{2BTfie), 
with 6ml = ml± — mlo, are made. In SU{2) xPT this contribntion appears at N^LO. 

The one-loop diagrams in Fig. give also the isospin-invariant nncleon-pNG conplings 
defined in Eq. (36). For these conplings we work in the isospin limit and find 


Vasgl = 


o'” 

yONEK 


Bm^d 


mx — QDF + 9F^) — 9{D -|- F)‘^m^ + {D — 


2 C2 

^IT 

Bm^O 


{f'{mK,^) - 8/'(m^, A)) 


+ 18DF — 15F^ -|- mKm.„ + ml 


A8tiF^ 


mK + m-n 
2 


{D - F){D + 3F) ml + m^^mx + m^ 


487rF2 


mx + mri 




2887r2F2 


A) -/(m^, A) f{mrj,A) - f{mx,A) 

ru o o ~r o o 


m 


K 


mt 


m, 


yONAK ~ Bm^d 


3{D‘^ — QDF — 3F^) m\- + m^^mx + m 


m 


K 


IGttFq 


mx + ^^71 

/72 


(59) 


(60) 


{D - F){D+ 3F)ml + m^mx+ mj^ f{mx, A) - f{m^, A) 


WtiF^ 


mx + m^i 


+ 


87r2F2 


m 


K 


mt 


(61) 


A. Testing the relations at NLO 


Armed with the NLO expressions for the baryon masses and the^ nncleon-pNG conplings, 


we investigate the relations fonnd in Section III We start with g^, which is of the largest 


phenomenological interest. We repeat the relations we want to test 

mj 


go = 6m 


N- 


me 


= (Am^ — Ams)- 


2m* 


-e 


(62) 

ms — m 

As these relations hold at LO, it is snfficient to test the relation for the NLO corrections 


themselves. A comparison of Eq. (45) and Eq. (56) shows that it is possible to write 


= 25 


6Wm^ 


N 




1 

+ 3 


73 


-1 


5(1) 


m 


7] —77 

'N 


ml — ml 


m^d 


(63) 


where 6^^'^m^ {6^^">m\ ’^) denotes the pieces of Eq. (45) induced by the kaon mass splitting {g 


- TT mixing). Using the LO expression for the meson masses and mixing angle this simplihes 
into 

= 5 PW 7 V . (64) 


me 
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Thus, NLO corrections, both with octet and decuplet intermediate states, conserve the relation 
between go and 6m n ■ 

Next we consider the second equality in Eq. (62). Using Eqs. (53) and (54), and expanding 


for simplicity the decuplet contributions in the limit A —0, the NLO corrections to go can 
be expressed as 


go = ( — A*-^^ms) + 


3(L)2 _ Qjgp _ 3^2) ^ 

2887rF2 


{mx - mT,Y{mK + 2m^) 


2m* 


m.. 


_-e . 

m 
(65) 

So in addition to a term proportional to Adlrnc — A*^^Ws, there is a second term that violates 
the relation. This second term vanishes in the SU{3) limit, m^ = m, and is non-analytic in 
the quark masses. The severity of the breaking is best illustrated by plugging in numerical 
values. Up to NLO it is possible to write 


9o 


Am= — Ams 


= 1 + 


3{D‘^ — 6DF — 3F^) + (mx — mnY{mK + 2m.n-) 


= (1-0.7+ 0.2) 


2887rF2 

2m* 
ms — m 


Ams — Ams 


2m* 


mg — m 


-9 


-0 , 


( 66 ) 


where the second and third contributions in the second line are the octet and decuplet correc¬ 
tions, respectively. We used the observed value of the S — S mass splitting, Am^ — Am^ = 124 
MeV [2|, and Fq = iA, the difference being higher order. We see that the tree-level relation is 
violated by a 50% correction and thus it is unsuitable for a precise determination of go. 


Next we look at the ^ couplings in Eq. (36). Using the NLO results for octet masses, we 


conclude that the following relations survive NLO corrections: 




o'" - 

yONSK ~ 


a 


( 1 ) 

Ns 


m.o 


i I . 


^9onak = - - A^^Wat) 


2m* 6* 


m ’ 


m^ 

6m* 0 
mg — m 


(67) 

( 68 ) 

(69) 


The remaining LO relations in Eqs. (40), (41) and (42) are violated. We observe that NLO 


corrections do not spoil the relations if the baryons that enter the ^ vertices are the same as 
those appearing in the mass combinations, while relations to masses of baryons that are not 
involved in the ^ vertices are violated. 

Finally we discuss g^. At tree level, gi is closely related to gor/, 9i^^ = 9o^,j(/> as can be seen 


from Eqs. (35) and (40) . This can be understood because at this order gi is induced by the 


emission of a r; meson by the nucleon and the consequent mixing of the g with a neutral pion. 
As was the case for g^]^, the hrst two relations of Eq. (35) are violated at NLO, however, in 


this case also the third relation is violated. At this order, gi receives a contribution identical 
to gQ^^, but in addition there are contributions from the kaon mass splitting and from g - 
7T mixing in the internal pion and g propagators. Neglecting in this discussion the decuplet 
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correction, we find 






5D^ — 6DF + 9F‘^ 
327iF^ 


{rriKO — mK+) 


1 (£>2 + qdF - 3F^){mr, - + 6{D^ + 3F‘^){ml - ml) 


247rF2 


mri + mj, 


Vs 


(70) 


The first piece, proportional to gorj, respects the relation to the nncleon sigma term. On the 
other hand, we were not able to hnd any nsefnl relation respected by the remaining piece of 


Eq. (70). As nnmerically the violation of the relation is of similar size as the g^^(j) part, the 


tree-level relation is of little nse. 

To assess the importance of corrections to the tree level valne of gi, we evalnate Eqs. (31) 
and (57) nsing bo = 0.068 GeV“^ and bp = —0.209 GeV“^ (these valnes are discussed in 
Sect. 


VIB). We hnd 


9i 

2K 


= (0.85 + 1.0 + 0.85) ■ 10-^0 


(71) 


The hrst number is the LO contribution. The second and third numbers are the octet and 


decuplet contributions to the NLO corrections, Eq. (57). We see that NLO corrections are 


large, as big as the leading term, and the inclusion of the decuplet makes them even larger. 
Since there is no surviving relation between gi and baryon masses or sigma terms, and the 
xPT corrections show no sign of convergence, we conclude that in SU{3) yPT there is no safe 
way to extract gi from the baryon spectrum. 


V. OCTET BARYON MASSES AND T' COUPLINGS AT N^LO 


In Section |IV| we have seen that NLO corrections affect the octet baryon mass splittings 
and nucleon ^ couplings in such a way that most LO tree-level relations beween ^ couplings 
and baryon masses are violated. The exceptions are the relations between g^ and the nucleon 
mass splitting Smjsi and similar relations for the couplings goNiiK and goNAK to Am^ — Am^ 
and Attia — Am^f. Furthermore, the link between ^o?? and the nucleon sigma term also 
survives NLO corrections. Of these couplings, go has the largest phenomenological impact as 
it contributes to the nucleon electric dipole form factor at LO an d give s ris e to t he dominant 
piece of the i9-induced ^ nucleon-nucleon potential. In Sections VA and VB we therefore 


focus on and go- We discuss gi in Section |V G[ and the nucleon couplings involving rj 
and kaons in Section IV Dl 

We show that again all loop corrections affect bm n and go in the same way, except for finite 
contributions that are quadratic in the isospin-breaking parameter e: and thus numerically 
suppressed. All LEGs from that are needed to absorb divergences in the loops appear 
in the same way in go and 6m]\f. The relation between go and is broken by additional, 
hnite contributions to go from operators in which do not contribute to dm^. These 
contributions, however, scale as B'^fh? and not as B'^msfh. Even though the values of these 
LEGs are currently not known, they should not affect the go — 6m n relation in a signihcant 
way. Finally, S'17(3)-breaking effects induce subleading pion and g tadpoles, which contribute 
to go but not to 6mpf, thus violating the relation. We estimate these violations and show that 
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(a) 


(b) 


(c) 


(d) 






(g) 



FIG. 2: 0{q^) corrections to the baryon octet mass splittings. Circled dotted vertices denote SU{3) 
invariant couplings of Eqs. (22) and (23). Crosses denote insertions of the octet and decuplet baryon 
mass terms of Eq. (20). Other notation as in Fig. 


our results allow for a precise extraction of go from available lattice evaluations of the strong 
part of the nucleon mass splitting. 

In what follows we calculate N^LO correction to the nucleon mass splitting (the nucleon 
average mass and the masses and mass splittings of the other octet baryons are given in 
Appendices]^ andand go including corrections due to the baryon decuplet. We keep terms 
linear in the quark mass difference, e, neglecting contributions. We comment on these 

corrections briefly at the end of the section. 


A. Corrections to the nucleon mass splitting 


N^LO corrections to the nucleon and Delta mass splittings in SU{2) xPT were considered 


in Ref. 


In SU{3) xPT N^LO octet masses and mass splittings were considered 


m 


Ref. [75], in the infrared regularization scheme. Here we repeat the calculation in the heavy 
baryon formalism, and include corrections from the baryon decuplet. At 0{q^), baryon masses 
receive corrections from loops involving vertices in the Lagrangians and in Eqs. 

(20), (22), and (23), and from tree-level insertions from the Lagrangian in Eq. (24). The 
relevant loop diagrams are shown in Fig. Diagramsj^a) - (f) show the contributions of octet 
intermediate states. Diagram j^a) includes the correction to the propagator. It contributes 
to the mass splittings in two ways, through the kaon mass splitting or g - n mixing, and 
through the on-shell relation, which relates v ■ p to the mass of the external baryon in the 
diagram. Diagrams [^b,c) contain recoil corrections to the axial couplings D and F, and are 
proportional to l/ms- Diagram j^d) has a piece given by a recoil correction, and a piece 
proportional to the LECs bi, 62 , and 63. The LEG 63 does not contribute to mass splittings, 
but only to isospin-averaged masses. Diagrams [^e,f) have an insertion of the SU{3) breaking 
couplings bn, bp, and bo- 

The diagrams in Fig. [^are UV divergent and the divergences are absorbed by the coun¬ 
terterms in Of the operators dehned in Eq. (24), only di, d 2 , d^, d^, and do are relevant 
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for mass splittings. Furthermore, di satisfy the relation in Eq. (26), implying that there are 
only four independent counterterms. We write 

Sm^j^ = {4B)‘^{2me){rhdi + msd 2 ) , 

Sm^ = {4:BY{2The){md3 + irisdi) , 

= —{4:BY{2me) {ffi{d^ + di) + 


6m^ 


d^ + (^2 


(72) 


where we neglected terms of order and di are dehned as 

di = —2((ii d2 d;i d^ dg) , d2 = 2di — 2d3 — dg — dg 

dg = 2(di — d2 dg — dg -|- dg) , d4 = — 2d\ 2d'^ — dg -|- dg 


(73) 


Because the counterterms satisfy the Coleman-Glashow relation, the divergences of the dia¬ 
grams in Fi g. must do so as well. We have explicitly checked that this holds and that, at 
N^LO, Eq. (26) is only violated by the finite term UTI 


6misf + 6ms + 5ms = 


^ 2hDDF) 




m., — m)mj^ log —^— 

'^K+ 


m 


2me m^ log 


m 


K 


mt 


B(Qbji) -|- 126 ^ -|- 7bc) 
367r2F2 


{{ms - m)f^ {mx, A) - me {f:^{mx, A) - 2/2(m^, A))) , (74) 


where the decuplet loop function f 2 {x,y) is dehned below in Eq. (81). 

We now present our results for the nucleon mass splitting, after having subtracted the UV 
divergences in the MS scheme. In order to facilitate the comparison to loop corrections to Qq, 
we split the N^LO corrections to 6mx in three contributions. We start from diagrams [^b,c,d) 


and the piece of[^a) proportional to 


m 


-1 

B 


= {bi + b 2 + 63 ) g 2 F 2 ( '^KO - ^K+ + ^KO log -m^+ log 

STT \ m^o 

1 




m 


K+ 


(36i + b 2 - 63 ) 






4 , 4i A 4, A 

m„ m^ log — - m log — 


- QDF - 3F2 


m 


4 1 „ 

RO 


log 
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mt 


m, 
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K+ 


log 




RO 


m 


K+ 


{D-3F){D + F) 0 / 4 , 

-—-^ log — - m„ 


mt 


.log-^ 

m2 


167r2Fo2mB ^/3 \ 

Here and in the following we omit the superscript (2). Diagrams [^a,e) give 

6mA'' = —8Bme{bx + bF) 


(75) 


167r2F2 


3D^ + 2DF + 3F2 
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-4:Bme{bi:) — bp) 


{D + FYui: 

A 




1 + 3 log — 


mt 


m^ 1 -h 3 log 


m 


K 


ms — m 
3Qffie 


{D + 3Ffg2{mxo,mx+) 
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127r2Fo2 


{D - F)m\ ( 1 3 log 
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ms — m 
6fhe 


{D - 3F)g2{mxo,mx+) 


(76) 
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where pion loops are proportional to the combination hr, + hF-, which determines 5m n at tree 
level, while kaon loops also give contributions proportional to hr, — bp- Only pion and kaon 
loops contribute to In the case of the rj and oi rj - tt mixing, the contributions of 

diagrams |^a) and [2|(e ) exactly cancel. 

Finally, diagrarnMf) gives 


5m^^ = —8Bme {bp, + bp) 


327r2F2 


m„ 




1 + log — + ^ 1 + log — 
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[mg — m I 
,2 


3 

2 2 I 2 1 2 1 

m^-m„ + m^ log — - log — 


-r, 


mt 


m. 


+m^ 1 + log 




m 


K 


m + ms , , 

+ 2m£ diimpo, mK+) 


mt 


The loop functions gi and g 2 are dehned as 

/^2 ^2 

gi{mKO,mK+) = m^o - m^+ + log -m^+ log , 

m^o rnj^+ 

g 2 {mK 0 ,mK+) = m^o - + 3m^o log - 3m^+ log 


(77) 


(78) 

(79) 


Eqs. (75), (76), and (77) reproduce the results of Ref. 

We then consider f^LO decuplet corrections to the baryon octet mass splittings. The 
relevant diagrams are depicted in Fig. [^g,h). We give here the results for Smpf, while the 
contributions to the S and S mass splittings can again be found in Appendix 

The UV poles of the diagrams involving decuplet intermediate states have the form 


/^''(m. A) = — (2A^ — m^) L . (80) 

The divergence proportional to A^ is absorbed by the 0{/S?) piece of the LECs bp and bp, 
while the divergence proportional to the quark mass is cancelled by the counterterms di and 
d 2 . After subtracting the UV poles, the decuplet contributions can be expressed in terms of 
the function 

,,2 _ A 

/ 2 (m,A) = —2A^ + (m^ — 2A^) log—- — AA\/m‘^ — arccot^ , (81) 

m2 y/fn2 _ ^2 

and we define 

fti.'rriK, A) = {f2{mK0, A) ± /2(mi^+, A)) . (82) 

Diagraming) induces corrections to the mass splittings that are proportional to the LECs b^, 
bp, and bp. 


5m^^^ = --^^^Bme{bp + bp) (8/2(m^, A) +/2+(mj^, A)) 

+ 247 r 2 F 2 ^ (("^^ + 2m)bo + m{bp + bp) + ms{bp - bp)) {mp, A) , (83) 
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where we neglected relativistic corrections to the decuplet propagator. Diagram [^h) contains 
corrections that are induced by mass splittings of the decuplet, and are proportional to the 


LEG be in Eq. (20) 


brn^N = 


heC^ 

'367r2F2 


B I 2me (l0/2(m,, A) + f}{mK,l2.)) - 


(84) 


Contributions proportional to 6 a in Eq- (20) can be absorbed by a redehnition of A and we 


did not explicitly consider them. If we neglect kaon loops, Eqs. (83) and (84) agree with the 
SU{2) calculation of Ref. [67] . 


B. N^LO corrections to 


go 


The N^LO corrections to with octet and decuplet intermediate states are shown in Fig.[^ 
In addition to these diagrams, we need the nucleon wave function renormalization Zn, and the 
pion wave function renormalization, Z-j^. For the calculation of go, it is sufficient to compute 
them in the isospin limit rriKO = 'mK+ = rriK and 0 = 0. Denoting SZ^ = Zn — 1 and 
SZt^ = Zjt — 1 , we have 


6Zn = 3{D + Ff 


mt 


647r2F2 


mt 


1 + 3 log -^ ) + (5D^ - QDF + 9F 


2^ 


967r2F2 


1 + 3 log 




m 


K 


+ {D-3F) 


2 




(85) 


5Z^ = 


ml 


2A7t^F^ 


1 + log At 


m 


K 


m‘i) 487r2Fg 


1 + log 




m 


K 


16R 

~w 


((rUs + 2m)L4 + mLs) . 

( 86 ) 


Finally, at this order one has to consider the correction arising from expressing Fq in the LO 
,2^ pion-nucleon coupling in terms of At N^LO, the relation between Fq and is 1121 


SF^ = l-f 

-CO 

2 

mi 


167r2F2 


1 + log - 
mi 


m 


K 


327r2F2 


1 + log ) - -f^{{ms + 2m)L^ + rfiLk) 


p2 

-CQ 


(87) 


We compute the diagrams in Fig. with on-shell baryons, and an incoming pion with 
energy v ■ q and zero three-momentum ^=0. The on-shell condition can be written as 


v-K = -[v-p + v-p) =-[mi + mf - 2mB) i -r- 

1 1 2771B 


V ■ q — {v ■ p' — V ■ p) = mf — nii , 


( 88 ) 

(89) 
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FIG. 3: N^LO contributions to ^ nucleon couplings. Squares denote ^ pNG and baryon-pNG 
couplings. Other notation is as in Figs. and . 


where ruf^i are the corrections to the masses of the baryons in the hnal and initial state, and 
rriB is the common mass of the octet. For the diagrams in Fig. [3, it is enongh to use the 
tree-level expression of the baryon masses, which are given in Eqs. (2^ and (27). 

The counterterms are determined from Eq. (24). Operators di - dg contribute in the 
same way to go and 6m^ and thus preserve the tree-level relation. Some of the remaining 
counterterms do spoil the go - 5m n relation. These corrections are discussed in Section fV B 1 
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We move on to the loop diagrams, which we discuss in some detail. The contribution of 
diagrams |^a,b,c,d) to % can be written as 


^0 


(a,b,c,d) r (a,b,c,d) r. 

0 = ora'^N ) 


me 


(90) 


where we applied the LO expressions of the meson masses and mixing angle. In this expression 


5mis given in Eq. (75). Both in Eq. (75) and in Eq. (90), only the pieces of diagrams 


l^a) andg a) proportional to 1 /ms are considered. 

When we combine the contribution of diagram [^a) proportional to v ■ K, that on-shell 
becomes proportional to the average nucleon mass, the nucleon wave function renormalization 
Z]y, and diagrams |^e,g), we obtain 


;;(a,e,g) 

VO 


-f gl)°^SZjv = —6 

me 


(a,e) 


(91) 


where is given in Eq. (30) and in Eq. (76). Contributions from the rj meson in the 

loop diagrams [^e,g) are cancelled by the sum of wave function renormalization and diagram 
[^a), in the same way as it happens for the nucleon mass splitting. Diagrams [^k,l) mutually 
cancel. 

Then, we consider diagrams ,j). These diagrams need to be combined with one- 

loop corrections to the pion wave function renormalization and to the decay constant Fq. 
Considering these effects, we hnd 






6F^ 1 = 5m® , 

me 


(92) 


ff) A -k 

where 5my is in Eq. (77). Notice that the contribution of the LECs L 4 and L 5 cancels 
between the pion wave function renormalization and SF^^. 

Finally we consider the decuplet corrections in Fig. [^o,p,q). An explicit calculation 
shows that diagrams |^o) and [^p,q) are in direct correspondence with Fig. |^g) and[^h), 
respectively, once the decuplet corrections to the nucleon wave function renormalization are 
included. Thus N^LO decuplet corrections do not spoil the go - 5m n relation. 

The results in this section show that all UV-divergent contributions to go up to N^LO can 
be expressed in terms of the strong part of the nucleon mass splitting. However, at this order, 
we hnd some hnite violations of the relation which we discuss now. 


1. I^LO violations of the relation between go and dm^- 


There are three types of N^LO corrections to go that cannot be written in terms of Sm^. The 
hrst of these corrections arises from additional counterterm contributions to go- As discussed 
in the last section, the counterterms di to dg in Eq. (24) conserve the relation, while dg - die 


can potentially spoil it. It can be seen that dg and di 2 do not contribute to ^ baryon-pNG 
couplings at this order, while dis and die only contribute to gi- The remaining correction can 
then be written as 

- {ABy de {2m) m^ 9 , (93) 


9o = 5m% — 


me 
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where is the combination 


— 2dio + 4(iii + 3(ii3 + 3d 


il4 


(94) 


The combination of LECs d^ is thus not related to mass splittings in the baryon spectrum, 
but could in principle be extracted from a precise analysis of nucleon-pion scattering. In 
practice, however, these LECs appear at too high order, and are not well constrained |78j . 
This additional counterterm is present in SU{2) yPT as well, where it also appears at N^LO 
Both in SU{2) and SU{3), the additional contribution to Qq scales as rh and not as 


in contrast to the terms in Eq. (72). Considering the good convergence of SU{2) yPT, we 
expect these corrections to be of the expected size, m^/A^, of the order of a few percent. 

The second type of contributions that violates the go - Sm^ relation appears due to addi¬ 
tional isospin violation. We have not calculated these contributions systematically, but give 
just one example. Diagramj^a) induces, besides the component proportional to v-K discussed 
above, a contribution proportional to the energy transfer v ■ q. Such a term gives rise to a 
correction to Qq that is quadratic in the quark mass splitting which is not matched by an 
analogous correction to hm^r. This correction is proportional to the tree-level pion-nucleon 
coupling, and we can write 


r-(c) 

_( 0 ) 

9o 


- 6DF - 3F 


2 m 


KO 


m 


K+ 


+ 2m\^m]^+ log 


K+ 

A 


2567r2FQ 


m'^KO - 


Plugging in the values of the axial couplings and kaon masses, Eq. (95) gives 


r-(c) 

-( 0 ) 

% 


~ -2 ■ 10 


-6 


(95) 


(96) 


which is completely negligible. We do not expect significant corrections from the remaining 
0{e^) contributions that we did not compute. 

Finally, the third type of violations arises from tadpole contributions, diagrams |^m,n). 
In the 0{q'^) meson Lagrangian, the operators L 7 and Ls generate pion and r] tadpoles. In 
addition, one has to consider one-loop diagrams with one external rj or no, and the ^ three- 
pNG vertex of Eq. (12). Together with the two-pNG vertices from the baryon mass terms, 
the tadpoles generate contributions to go and gi. Up to corrections of the contribution 

to go arises only from rj tadpoles and is proportional to We hnd 


S9i 




-( 0 ) 

9o 




—3m 


32B 

F2 

-^0 


(Lg 3Lj){ms -m) + 


32n^F^ 


2m\ log 


A 


m 


K 


2 1 A 
+ m log ^ 


2 1 A 

log — 
mi 


(97) 


where here we have adopted the subtraction scheme of Ref. m to dehne the renormalized 
coupling Lg, while L 7 is not renormalized. Using the values of L 7 and Lg discussed in Section 
II A[ we can estimate the tadpole corrections to be 




^9o 


-( 0 ) 

9o 


= (-0.5 ±1.3)-10 


-2 


(98) 


We thus expect the relation between go and 6m]\f to hold up to a few percent. 
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C. N^LO corrections to gi 


We have shown that already at NLO no relations between gi and baryon masses survive. 
However, in SU{2) yPT, where gi only appears at N^LO, it can be shown that part of gi 
can be related to the pion mass splitting induced by the quark mass difference [H]. At the 
same order, there appears an unknown direct contribution to gi which has been estimated to 
be small by use of resonance saturation techniques in Ref. 173 ]. In this section, we study this 
relation in SU{3) yPT and study the effects of the SU{3) LO and NLO contributions to gi 
that are missing in SU{2). 

The SU{2) relation of gi to the strong component of the pion mass splitting can be recovered 
by studying the tadpole diagrams. Fig. (ni,n), and, in particular, the contribution of the 
pion tadpole. In SU{3) yPT the pion tadpole receives contributions from the LECs L 7 and 
L^, and from one-loop diagrams with insertion of the^ three-pNG vertex in Eq. (12). When 
these contributions are combined, the coupling of the pion to the vacuum given by 


Aad = -2BmJ /tad FqTI^ 


(99) 


where the function /tad is 
48R 


/tad - ( 3^7 + LI) + 227 ^ 2^2 


2 1 

rrij^o log 


m 


V332n^Fl 


2m? 


K 


log 






2 1 F 

m^+ log — 


RO 


m 


K+ 




m 


+ 3m log — - 5m^ log — 


K 


mt 


mt 


( 100 ) 


It is possible to show that /tad is related to the pion mass splitting in the large limit, more 
precisely 

6m? 

lim /tad = - lim —^ , ( 101 ) 

ms—^00 ms—^co m^E 

where 6m\ denotes the component of the pion mass splitting induced by m^ — The 
expressions of the pion mass and mass splitting at one loop can be found in Refs. [l2l [79] . 
The pion tadpole induces a correction to gi, of the form 


= -8B(2b„ + 6 b + bF)ft^m,e , 


( 102 ) 


which, using Eqs. 
and Sm^ as 


(27) and (101), can be expressed in terms of the light quark sigma term 


r-(m,n) 

^^1 


m-—A^^^m]sf 

dm 


5ml 


mt 


(103) 


where we expanded also m*/m in the large limit. Eq. (103) is exactly what is found in 


SU{2) [Hj. Notice, however, that in SU{3) it is not possible to express gi in terms of the full 
pion mass splitting, but only of its large limit. 

To estimate the pion tadpole correction to gi, we use the extraction of the light quark sigma 
term of Ref. |80]. In this paper, the SU{3) yPT expressions of the nucleon mass and sigma 
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term are fitted to lattice data, and the LO contribution to the light quark nucleon sigma term 
is found to be 

... ( 104 ) 


m- 


dm 


= 65 ± 19 MeV 


Using Eq. (104), together with L-j and Lg discussed in Section II A, the pion tadpole contri 


bution to Qi is 


5g 


(m,n) 


2F. 


= -(1.9 ±3.7) ■ 10-^0 


(105) 


The central value is in reasonable agreement with the estimate of Ref. [22], while the larger 
errors stem mainly from the uncertainties on Lj and Lg, and the partial cancellation between 
Li and Lg in the combination SLy ± Lg. 

Eq. (103) is only a subset of the N^LO corrections to gi. In addition one should consider 
the loop diagrams in Fig. the r] tadpole, contributions from corrections to p - vr mixing, 
and N^LO counterterms. We have not computed these corrections, since, as we discuss now, 
they are not particularly instructive. 

Comparing Eqs. (105) and 0. we see that the contribution of the pion tadpole, which 
is formally N^LO, is comparable with the LO and NLO pieces of gi. This is not surprising, 
since the LO and NLO terms vanish in the limit —)■ oo, and thus are suppressed by powers 
of fh/rris. On the other hand, at N^LO gi starts to receive contributions that are hnite in the 
rris —>■ oo limit, pion tadpoles being one such example. Eqs. (71) and (105) show that for 
gi the suppression due to inverse powers of rris or is similar. An analogous observation 
was made for which is also determined at LO by the g - tt mixing angle m- Thus, we 
conclude that the SU (3) yPT power counting does not provide a good organizational principle 
for gi. Currently the best possible estimates are those based on SU{2) [22], which, however, 
are also affected by large uncertainties, leaving the determination of this important coupling 
in an unsatisfactory status. For clarity we repeat the SU (2) estimate of Ref. [22] here 


-SU{2) 

9i 

2F^ 


= -(3.4 ±1.5) ■ 10-^0 


(106) 


which is partially based on a resonance saturation estimate of an unknown N^LO LEC 
and on the tree-level relation Sm^. = /{rn\- — rn\). The error in Eq. (106) is perhaps 

slightly underestimated. Including the uncertainty on the pion mass splitting, either through 
the estimate of higher-order corrections [12] , or using the error on the extraction of Ref. 


6ml = (87±55) MeV^, would raise the 45% uncertainty in Eq. (106) to about 70%. We stress 


that even with this large uncertainty, the SU (2) determination of gi is incompatible with the 
LO SU{3) tree-level estimate in Eq. 0- Finally, we note that Eq. ( |106[ ) does not include yet 
a higher-order contribution arising from the ^ three-pion vertex^. This contribution enhances 


the estimate in Eq. (106) by roughly 50% [22l [28] making the tree-level estimate of gi even 
less reliable. 

The large uncertainty on gi dominates the uncertainty of light-nuclear EDMs when ex¬ 
pressed in terms of 9 directly [33]. It also affects the EDMs of diamagnetic atoms, but there 


^ In SU{3) yPT this contribution already appears at NLO and is identified with the piece proportional to 
{D + F)^m^ in Eq. ([^). 
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the nuclear uncertainty associated with the complicated nuclear many-body problem is still 
larger. Nevertheless more theoretical work on the size of gi could signihcantly increase the 
precision of EDM analyses. 


D. Other ^ nucleon couplings 


The relations of goNJ^K, and goNAK to the nucleon sigma term, the nucleon-S and 
nucleon-A mass splittings can be checked in a way analogous to what done in Sections V A 
and |VB| for qq and 6mN- The expression for ArriN, Am^, and Am\ are given in Appendixj^ In 
addition, one needs the expressions of the baryon and meson wave function renormalizations, 
and the expression of and as function of Fq, which we also give in Appendix]^ After 
all the ingredients are put together, it is possible to show that the loop contributions to the 
X couplings satisfy 


y^gov = + gi°rl + 


SF^ 1 = 
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( 2 ) 

Ns 
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( 2 )' 

Nl 


m. 


2m 


4m* 0 


( 107 ) 


gO NEK — go NEK + go NEK ( 2 SFk 


=(loop) 
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= — (A‘'^Ws — A'-^Wat) 


2m^9 
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'^goNAK — 


-(loop) . _(0) 

gONAK + gONAK 


- {6Zn + 6Z^ + 5Zk) 


SF, 


K 


= — (A'^^Wa — A*^^^mAr) 


6m^6 


m, — m 


(108) 


(109) 


Eqs. (107), (108) and (109) include all the loop corrections, with the exception of the v ■ q 
contribution to diagram |^a). In the case oi go^, v ■ q = 0, and this contribution vanishes. 
For go nek and goNAK, this contribution violates the tree-level relations, as in the case of go. 
However, for go the breaking is proportional to while for couplings involving the nucleon 
and the S or A baryons it goes as (m^, — m)^, and thus is potentially larger. We hnd 


Sgo NEK,v-q 
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5127r2F2 


(F)2 - 18DF + 9F2 
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+9(D - Ff 


2m\ml log ^ 
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5go NAK,vq 
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gONAK 


bUir^F^ 
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( 111 ) 
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Therefore, also for goNT,K and goNAK, these violations are only a few percent. 

The other contributions that violate the tree-level relations are the tadpole diagrams, 
[^m,n), and the counterterms dg - diQ. For the nucleon-S and nucleon-A couplings, the 
tadpole contribution is proportional to the tree level, and we can write 


-( 0 ) 

doNEK 


NAK 
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9o NAK 
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( 112 ) 


where, again, we have used the subtraction scheme of Ref. |12]. Using the values of L-j and Lg 
discussed in Section II A , we find that the tadpole corrections amount to no more than 10%. 
In the case of the g, the tadpole corrections are not proportional to the tree level 
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(113) 


In this case, the estimate of the tadpole corrections is affected by larger uncertainties. Nonethe¬ 
less, using the tree-level values of bo and bp, bp = 0.068 GeV“^ and bp = —0.209 GeV“^, and 
expressing 2bo + bp + bp in terms of the LO contribution to the light quark nucleon sigma 
term, Eq. (104), we get that the tadpole corrections to ^or? come in at 30%. 


Finally, the operators in the Lagrangian can also violate the tree-level relations. The 
counterterms di - ds, which are needed to absorb the divergences in the baryon mass, respect 
the relations, as one expects. The operators dg - dn only contribute to the T couplings, and 
we hnd 




0 Tj 


m^6 


= -(4R)^(6dio - 4dii 9di3 - Odu) -{2ms + m) 

o 

-|-(4i?)^8 ( —dll + 2(ii4 -|- 3(ii5 -|- dig j {mg — fh,). 




0NY.K 


m^6 


b9, 


ONAK 


m^6 


= (4R)^(2dio - 4dii 3di3 - 3di4) (m^ -1- m), 


— (4i?)^(6dio T 4dii -|- 9di3 -|- 3di4) {mg -h m). 


(114) 

(115) 

(116) 


We cannot give a precise estimate of the corrections in Eqs. (114)-(116), due to the ignorance 
of the LEGs dio, dn, and di3 - dig. However, we see that the corrections to gori scale, at 
most, as m^/A^, and to goNj:K and (Jonak as mj^/Al- Taking A^ = dTrU^-, the counterterm 
corrections should be around 15% - 20%, of similar size as the tadpole contribution. 
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VI. DISCUSSION 


A. Best estimates of the 'X couplings 


In this section we give the best estimates of the X couplings by using the relations in Eqs. 


(33) and (67)-(69) that do not suffer from large SU{3) yPT corrections. 

The coupling go is related to the strong contribution the nucleon mass splitting for which 
there are now several lattice QCD calculations [271 1321 - 186] . The first three calculations [821 - 
IM] were performed with only a single lattice spacing and with pion masses m.,r ^ 250 MeV. 
According to the FLAG Lattice Averaging Group standards HU, these results would not be 
included in averages of lattice QGD predictions. The next calculation [85] was performed with 
> 283 MeV and four lattice spacings. According to the FLAG criterion, these results would 


be included in an average, but perhaps do not have all the systematic under complete control. 
The final two calculations EZllSS], although performed by the same group, are independent 
from each other, are performed with multiple lattice spacings and with pion masses at or near 
their physical value. Both calculations include effects on the splitting from QFD. In the hrst 
case, QFD was not included in the sea quarks, while in the second, the entire calculation 
included the effects of QFD. The strong contribution to the nucleon mass splitting from both 
of these calculations would receive the “green star” from FLAG. 

There are only a small number of results that pass the FLAG criterion. For a hnal deter¬ 
mination of the rud — rriu contribution to 6m^ from lattice QGD, one must have more results. 
In this case, the exclusion of the hrst three results mostly stems from the use of a single lattice 
spacing. In the mass splitting, the leading discretization effects exactly cancel, since the lattice 
regulator used in those works respects havor symmetry. We therefore chose to include all the 
results [271IH2HH6] to construct a lattice average. However, we assign a weight penalty to the 
hrst three calculations. This follows the averaging scheme in Refs. [261187] : the weight factor 
is chosen to be 


Wi = 



(117) 


where at are the given statistical and systematic uncertainties in the given lattice QGD cal¬ 
culation combined in quadrature and i/i = 1 for the hrst three calculations [H2HH1], l/i = 2 for 
Ref. [85], and yi = 3 for the most recent two calculations [271|86|, with these weights chosen 
somewhat arbitrarily. This weighted average yields 


6mN = 2.49 ± 0.17 MeV 


(118) 


Substituting this in Eq. (33) together with e, Eq. (|^, we obtain 

% 


2K 


= (15.5 ± 2.0 ± 1.6) ■ 10"'^ 0 


(119) 


The hrst uncertainty comes from combining those on dmjq and e in quadrature, while the 
second is a conservative estimate of the theoretical error associated to the N^LO corrections 


discussed in Section VB This estimate agrees with the recent determination in Ref. 


based on SU{2) yPT. Because the relation to the nucleon mass splitting is preserved, the 
only diherence compared to SU{2) is due to the m/mg correction in Eq. which is tiny. 
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due to inclusion of 


In addition, the error in Eq. (119) is slightly larger than in Ref. 
higher-order chiral corrections. 

The fact that the relations of goNiiK and Qonak to, respectively, mj^ — rriN and rriA — mN are 
only violated by finite N^LO corrections, allows for a reliable estimate of these couplings. In 
this case, the electromagnetic contribution to the isospin averaged masses is relatively small, 
and we can use the experimentally observed baryon masses, Am-s = 1193 MeV, Amjy = 1116 
MeV, and Atun = 939 MeV, with negligible experimental uncertainties. We use the ratio 
ms/fh in Eq. (Q to obtain 


go NEK 
2Fk 
go NkK 
2Fk 


= - (36±l±ll)-lO-'^0 
= - (44 ± 1 ± 13) ■ 10“^ e 


( 120 ) 

( 121 ) 


The first error is given by the errors on ms/m and £, while the second estimates the effects of 
the finite terms that break the relations. In this case the breaking scales as m \/which we 


estimate at the 30% level, see the discussion in Section VD 

Finally, gon is expressed in terms of the nucleon sigma terms. For the light quark sigma 
term, the most precise value is determined from low-energy kN scattering with the most recent 
determination from Ref. [88], 

cr^vi = 59.1 ±3.5 MeV . (122) 

This number is consistent with earlier yPT analyses [HnHSI]- This quantity can also be de¬ 
termined with lattice QCD. However, there is significantly larger uncertainty from the lattice 
determination arising from a few systematic issues. The primary means to determine this 


quantity is invoking the Feynman-Heilman theorem, Eq. (43), with a large spread of results, 
see Ref. [92] for a recent review. There is a surprising “phenomenological” pion mass depen¬ 
dence of the nucleon mass found in lattice QCD calculations [93l[9l| yielding rriN ~ 800 MeV 
RttItt over a wide range of pion masses, including the physical point. This in turns provides 
an estimate of = 67 ± 5 MeV. The best lattice QCD calculation, which also would receive 
a “green star” from FLAG was performed with pion masses as light as ~ 190 MeV [HD] 
with the result 

= 39+^® MeV . (123) 

It will likely be years before lattice QCD results can compete with the dispersive ttN scattering 
determination of Ref. 


in Eq. (122). 


For the scalar strange content of the nucleon, there is no close second to the lattice QCD 
determination, although the results are not yet mature. Ref. ETI compared all recent lattice 
QCD calculations of aNs [ZS] |80l I95lll01j (where only Ref. [80] evaluated all systematics) and 
found a systematically low value, as compared with prior estimates from SU{3) baryon yPT. 
An average value was determined 


o'ns = 40 ± 10 MeV 


Thus our estimate for the coupling using inputs from Eqs. (122) and (124) is 

90r] 


2F. 


= (115 ±8 ±35) ■ 1O"®0 


(124) 


(125) 


where the hrst error is from the uncertainty in the sigma terms and the second from the N^LO 
corrections discussed in Section IVDI 
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FIG. 4: An example of an NLO two-loop correction to baryon EDMs involving the 'J’ three-pNG 
vertex. Several other diagrams with different topologies appear at the same order. 

B. A comment on baryon EDMs 


The 0-induced EDMs of the baryon octet in the framework of three-flavored yPT have 
been studied in great detail in Refs. |2S1 SSI 11021 llU3j (for a calculation of the nucleon EDM 
arising from the CKM phase, see Ref. jlU4j ). In these works U{3)l x U{3)r yPT is applied 
to calculate the EDMs of the whole baryon octet up to NLO in the chiral power counting. 
Compared to the SU{3) x SU{3) yPT framework used in this work, the main difference arises 
from the dynamical inclusion of the r]' meson. However, as shown in Ref. [102] . these effects 
can be absorbed in a redefinition of a counterterm contributing to the EDMs of the charged 
baryons p, and 

In the case of CP violation induced by the 6 term, the first contributions to baryon EDMs 
arise at 0{q‘^/A^). At this order, the baryon EDMs depend on only two combinations of coun¬ 
terterms in addition to one-loop diagrams involving the ^ nucleon-pNG vertices go, goNUKi 
and goNAK, and analogous couplings of the other octet baryons to pions and kaons. At NLO, 
no new counterterms appear but several additional one-loop diagrams contribute with no ad¬ 
ditional LECs. At NLO, one finds the first contributions from neutral mesons, and thus from 
gor^ and gi. Corrections induced by gi, being proportional to the/mg, are small, and have been 
neglected in Refs. |23l 1102] . In SU{2) yPT, gi only contributes at N^LO |2l]. The baryon- 
pNG vertices are related to and bp using the tree-level relations in Eqs. (30), (37)-(39) 
and values for {bp, bp} = {0.068, —0.209}GeV“^ are obtained by fitting the baryon masses 


to the tree-level expressions in Eq. (27). The two unknown counterterms are fitted to lattice 


data of the neutron and proton EDMs and the EDMs of the other baryons are predicted. It is 
found that NLO contributions are significant, in particular for nonphysical large pion masses 
that are typically used in lattice evaluations of the nucleon EDM. 

One class of diagrams has not been considered in these works. This class consists of 
diagrams involving the ^ three-pNG vertices. Despite being two-loop these diagrams are of 
order and thus contribute to baryon EDMs at NLO^, and can be potentially large. 


^ These diagrams are NLO because the three-pNG vertices appear at lower order than the tree-level^ baryon- 
pNG vertices used in the LO one-loop diagrams. Although the extra loop comes with a suppression of 
(g/Ajj.)^, the relative size of the three-pNG vertex brings in a factor A-^/q making the two-loop diagrams 
genuine NLO corrections. 
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An example of such a two-loop diagram is shown in Fig. A two-loop calculation is beyond 
the scope of this work, but we note here that some of the two-loop diagrams involve a one-loop 
subpiece with the topology of Fig. [^c). Thus at least part of these corrections are taken into 
account by including one-loop corrections to go, Qq^, goNUK-, and goNAK, that is, they are 
automatically taken into account if the baryon-pNG couplings are obtained from the relation 
to baryon masses that survive higher-order corrections. Of course, a better assessment of 
the strangeness contribution to the nucleon EDM requires the full two-loop calculation. In 
two-flavor xPT ^ three-pion couplings arise only at subleading order, implying that two-loop 
diagrams as the one depicted in Fig. j^only contribute at N^LO and can thus be neglected. 
In Table I a comparison is given between the tree-level predictions in Eqs. ([^, (|^-([M|) 


using {bo, bp} = {0.068, —0.209}GeV ^ and the values obtained in Section VIA The two 


predictions agree well for the couplings to kaons. goNT.K and goNAK are only 10% smaller, 
and the predictions agree within errors. The coupling to pions is more affected, being roughly 
40% smaller, with a smaller uncertainty (this point was already made in Ref. [73] based on LO 
SU{3) yPT arguments). This can be easily understood, since using the values of bp and bp 
obtained by fits to the tree-level isospin-averaged octet masses is equivalent to use the relation 
of go to Arri’E — Ams, rather than the robust relation to 6m n. Finally, using the nucleon 
sigma terms rather than the tree-level prediction in terms of bp and bp leads to a considerably 
larger gor^- This is a reflection of the poor convergence of the SU{3) expansion for the nucleon 
sigma term. However, our best estimate is affected by a relatively large error due to unknown 
LEGs that enter at N^LO. 

We conclude that 517(3) corrections to the baryon-pNG couplings moderately alter the 
tree-level predictions. We do not expect that the 40% shift in go signihcantly affects the 
nucleon EDM extractions performed in Ref. [10211103] , as the lattice data are not yet sensitive 
to non-zero values of ^ nucleon-pNG couplings [25]. (Notice, however, that the analysis of 
Ref. [25] did not include the more recent results of Ref. [H].) For future extractions based 
on more precise lattice data with pion masses closer to the physical point, we recommend the 
values (and uncertainties^) of the ^ couplings given in the third column of Table 

Finally we discuss decuplet corrections to the nucleon EDMs, which have not been cal¬ 
culated in the literature. The leading^ Lagrangian induced by 6, Eq. (20), contains non- 
derivative^ decuplet-pNG and octet-pNG couplings, but no;2^ decuplet-octet-pNG couplings. 
Such couplings require, in order to conserve angular momentum, at least two derivatives, and 
thus are suppressed in the xPT power counting. The lack of a^ nucleon-decuplet-pNG ver¬ 
tex at the order we work implies that there are no LO one-loop contributions to the nucleon 
EDM. In three-flavor xPT; two-loop contributions appear at NLO due to diagrams with sim¬ 
ilar topology as Fig. but with an internal decuplet propagator. However, as argued above, 
these corrections are partially taken into account if the^ baryon-pNG couplings are inferred 
from the protected relations. We therefore do not expect decuplet corrections to play an 
important role in the study of baryon EDMs. We note that contributions from the decuplet 
EDMs to octet EDMs only appear at N^LO. 


^ A part of the uncertainties on the couplings arise from N^LO corrections and thus are formally higher order 
than considered in the baryon EDM calculation. They were estimated in Refs. [10211103j by varying the 
renormalization scale appearing in the chiral logarithms. 
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Tree-level values 
[x 10-3 0] 

Values obtained here 
[x 10-3 0] 

5o/(2FO 

26 

15.5 ±2.5 

90r)/ (20;) 

56 

115 ±37 

go nek/{2Fk) 

-41 

-36± 11 

go nak/ (2Er-) 

-48 

-44 ± 13 


TABLE I: Comparison between tree-level predictions for^ nucleon-pNG couplings using {b^, bp} = 
{0.068, —0.209} GeV“^ (see Refs. [23lll02] i. and the predictions from Section VIA All values are in 
units of 1O“^0. 


1. The nucleon Schiff moment 


An alternative way of extracting the^ nucleon-pNG couplings is to study the momentum 
dependence of the nucleon electric dipole form factor (EDFF). While the FDM gets both 
long-range and short-range contributions, the momentum dependence of the FDFF at NLO 
is hnite, and determined purely by loop diagrams. Short range effects enter only at N^LO. 
The full momentum dependence of the FDFF is given in Refs. [23l 1102] . with the omission 
of two-loop diagrams as the one depicted in Fig. Denoting the nucleon (proton) EDFF as 
, the Schiff moment is defined as 


S — — 

^n,p 


dFn,p{q^) 


g2=o 


and is given at NLO by 
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(127) 


(128) 


where we expressed the results in Refs. |23l 1102] in terms of go, go nek, and goNAK, and of the 
nucleon-S and nucleon-A mass splittings, and we used D + F = qa- Considering only pion 
loops, these results agree with the SU{2) calculations of Refs. [211 110511106] . In SU{3) yPT, 
the neutron and proton Schiff moments have two components. The contribution of pion loops 
is isovector and its typical scale is determined by the pion mass. In addition, both neutron 
and proton receive contributions from kaon loops, which vary on the scale of the kaon mass. 
The neutron only receives contributions from loops involving the S baryon, while the proton 
from both S and A intermediate states. One can immediately see that the kaon contributions 
receive large NLO corrections, indeed larger than the LO, being TimK/rriN ~ 2. However, the 
neutron and proton Schiff moments receive their largest contributions from pion loops, so that 
in this case the poor convergence is not a big issue quantitatively. 
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FIG. 5: Contribution of goNY:K and goNAK to the nucleon-nucleon^ potential. Single lines denote 
nucleon external states. Double lines denote A or S baryon propagators. Only one possible ordering 
per topology is shown. 


Using the values of the T couplings in Table [H we hnd 

Sn + Sp 


So = 


= - (0.5 ± 0.7) ■ efni'" , 


Si = = ((7.6 ± 1.4) - (2.4 ± 0.6)) ■ efm® , 


(129) 


where the isoscalar Schiff moment is given purely by kaon loops. For Si we have listed 
separately the contributions of pion and kaon loops. The uncertainties only include the errors 
from the ^ couplings, while we do not give an estimate of the theoretical error from higher- 
order corrections. Since the isovector Schiff moment is very sensitive to go, using our best 


estimate in Eq. (119) results in nucleon Schiff moments that are smaller by roughly a factor 

It is interesting that even in SU{3) yPT the Schiff moment is 


2 than found in Ref. 

predominantly isovector, due to accidental cancellations between the pieces proportional to 
goNT,K and goNAK- 

While a measurement of the nucleon EDFF is not going to happen in the foreseeable future, 
the predictions (127), (128), and (129) can be compared to lattice evaluations of the EDFFs. In 


particular, lattice calculations performed at several values of light quark masses could provide 
enough information to disentangle the contributions from the pion-nucleon and kaon-nucleon 
^ couplings, which have different dependence on rh and rris- This would provide a method to 
check the values of the couplings. 


C. A few comments on the ^ nucleon-nucleon potential 


The EDMs of light nuclei and diamagnetic atoms obtain important contributions from 
the ^ nucleon-nucleon potential. In case of the QCD 6 term, this potential is expected to 
be dominated by one-pion-exchange diagrams. In particular, because gi is suppressed with 
respect to go, often only the latter is taken into account. However, we stress that for certain 
^ quantities such as the EDMs of the light nuclei ^H, ®Li, and ®Be |in7flin9j . the ^^^Ra Schiff 
moment |8], or the ^®^Ta magnetic quadrupole moment [HQ], the nuclear matrix element for 
gi is signihcantly larger than that for go- In these cases gi should be included in the analysis. 
The value of go extracted in Section [Vl A| can be immediately used in the existing calculations 
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of EDMs that use the lowest-order chiral T potential induced by the QCD Q term [72] . On 
the other hand, as discussed in Section V C[ at the moment chiral symmetry only allows 
a determination of and thus of the T isospin-breaking potential, with a relatively large 
uncertainty. 

Nuclear EDMs can obtain important contributions from 'X nucleon-nucleon contact inter¬ 
actions. In case of strong CP violation, such interactions appear at N^LO only and are thus 
expected to be small. The corresponding potential is of the form 1721 


UrI?) = -f (Cl + C2T^'> ■ T®) (ff'" - ffW) ■ q 


(130) 


where and denote the spin and isospin of nucleon i, and ^is the momentum transfer 
q = p — p' 1 with p {p') the center-of-mass momentum of the incoming (outgoing) nucleons. 
Ci ^2 are LECs that in SU{2) yPT scale as Ci ^2 = 0{9rri^/C 2 is needed to absorb the 
divergences in two-pion exchange diagrams (TPE), while Ci is not renormalized at this order 
The TPE contribution can be used to give a rough estimate of C 2 


1 ^ 21=0 


(b(4rf.)= 


9o9a 


~ 2 ■ lO"^0fm 


(131) 


and Cl is expected to be of similar size. Clearly, this estimate is not very precise and below 
we derive an independent estimate of the sizes of Ci^ 2 - 

For nuclear physics applications, the typical momentum transfer is smaller than the kaon 
and r] masses. We can estimate the size of Ci ^2 by calculating contributions from gQ^, goNUKi 
and goNAK to the potential, and expanding them in powers of \q\/mK,ri, where \q\ ~ 171 -,^ In 
principle the same can be done for decuplet corrections but these only appear at higher order 
because of the absence of a^ nucleon-decuplet-pNG vertex (see the discussion at the end of 
Sect. IVlBl . 

The coupling to p, gon, contributes to the potential at tree level, providing a hnite piece to 
Cl- We obtain 


Ci = 


9oi 


D-3F 1 


2F2 

Tj 


V3 


= -8 ■ lO"^0fm' 


mt 


(132) 


which is in reasonable agreement with Eq. (131) considering the uncertainty of that estimate. 


The couplings goN'ZK and Xnak contribute to the nucleon potential only at one loop 
and are thus formally suppressed by two powers in the chiral counting compared to the one- 
pion-exchange contribution. Nevertheless, because the X coupling to kaons are somewhat 
larger than g^, their contribution might be sizeable. The corresponding diagrams are shown 
in Fig. 1^ They are similar to the TPE diagrams studied in Ref. [72|, with the exception 
that only triangle and crossed diagrams are possible because box diagrams are forbidden by 
strangeness conservation. We computed the diagrams in the limit \q\ <C rriK, while keeping 
the mass difference between the nucleon and the strange baryons in the baryon propagator. 
The diagrams contribute both to Ci and to C 2 - After cancelling the UV poles, diagrams 
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involving the T nucleon-S coupling give 


CiT. — 
C 2 T, = 
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SF^iAwFK)^ 
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while those with the ^ nucleon-A coupling 
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The function v{mK, A) is defined as 
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(135) 


and As = Ams — AmAi, Aa = Atua — AmAr. The first diagram in Fig. |^can involve E and A 
intermediate states at the same time, and thus is a function of As and Aa- In this case, we find 
that the loop function is well approximated by evaluating v in Asa = (Ams+AmA— 2 AmAr)/ 2 , 
and ignoring contributions proportional to Ams — ArriA- 

We can estimate the contributions to the potentials by setting the scale g, = rriN = 939 
MeV, and using the values for the^ couplings in Eqs. (120) and (121) 


Fis = -0.5-lO-=^0fm^ 
CiA = l.O-lO-^0fm^ 


C' 2 s= 0.1 ■ lO-'^fm'' 

C2 a = -0.5- 10-^0 fm=^ 


(136) 

(137) 


Thus, the contributions of the nucleon-kaon couplings to the ^ potentials are below the NDA 
estimate of Ci ^2 in Ep- (131), and signihcantly smaller than the contribution of ^or?- 

The contributions of Ci _2 to the EDMs of ^He and were studied in Refs. 1211 [22]. 
Unfortunately, these operators are very sensitive to the choice of the T-conserving strong- 
interaction potential and results vary by about an order of magnitude. Here we use the results 
of Ref. [22] obtained with the N^LO chiral potential |111] . which gives the largest dependence 
on C'i^ 2 - The EDM of ^He is, ignoring uncertainties, given by 


dsHe = 0.9d„ - 0.03dp + ( -0.11;:^ - OAOCiF^ + 0 

2F„ 


'CoFi I efm 


(138) 


Focusing on the pieces proportional to go and Ci, and using the estimates in Eqs. (119) and 

dsHe - 0.9dn + 0.03dp = (-1.8 + 0.3) ■ 10“¥ efm , (139) 


(132), we hnd 
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where the hrst (second) number is the contribution of go ((5i). The short-range potential 
provides a 15% correction compared to the one-pion-exchange contribution^, well within the 
0(20%) uncertainties of the nuclear calculation of the latter (see Ref. [22] for a more detailed 
discussion). It seems safe to neglect the short-range contributions to, at least, light-nuclear 
EDMs in case of the QCD 6 term. This conclusion is in line with Refs. f2T[ 122] . 


VII. CONCLUSION 


In this work we have investigated higher-order SU{3) flavor-breaking corrections to the 
relations between 2^ meson-nucleon couplings and baryon masses in the framework of SU{3) 
xPT. For each isospin-invariant 2^ nucleon-meson coupling induced by the QCD 6 term, we 
have identified one relation to baryon mass splittings or sigma terms which is not spoiled by 
loop corrections up to next-to-next-to-leading order. The determination of these couplings 
from spectroscopy is therefore independent of the convergence issues of SU{3) baryon xPT 
through this order. In Sect. VIA we have used the conserved relations to derive precise values 
for 2^ couplings. We recommend these values and corresponding uncertainties in future lattice 
QCD extractions of the nucleon EDM and in studies of nuclear 2^ quantities such as EDMs, 
Schiff moments, and magnetic quadrupole moments. 

The most important of these relations is the one linking the pion-nucleon coupling 2o to the 
nucleon mass splitting induced by rud — rriu, the quark mass difference. We find that all loop 
corrections, with the exception of small terms quadratic in rrid — mu, affect go and in the 
same way, so that at N^LO we can express go in terms of Sm^, plus corrections that are not 
enhanced by chiral logarithms, and are not proportional to m^. We stress that, at this order, 
the effects of strangeness on go are completely buried in Sm^, and thus accounted for when 
using lattice calculations of with dynamical strange quarks. We conservatively estimate 
the impact of terms violating the relation between 2o uud 6 m]\f to be about 10%. Thus, 
available lattice calculations of allow to determine go with 15% uncertainty. Further 
reduction of the errors will require improvements of the lattice calculations of the nucleon 
mass splitting, but also the determination of the unknown LECs that enter at Eq. 

(93), which are not related to baryon masses. Using the tree-level relation to Ams — Ams, 


as often done in the literature, overestimates go by about 50%. It will be interesting to see if 
a direct extraction of go from the lattice, e.g. from the momentum dependence of the nucleon 


electric dipole form factor, will give a value compatible with Eq. (119). 

Similarly, the 2^ couplings of the nucleon to the rj meson, and the couplings involving 
kaons, go nek and goNKK, are determined, respectively, by the nucleon sigma term, and by 
the mass differences of the nucleon and E and A baryons. In this case, the importance of 
the terms breaking the relation is larger, since they scale as m|^/A^ rather than m^/A^, and 
the corrections are estimated to be 30%. These couplings contribute to the nucleon EDM at 


LO, but, as discussed in Section VIC they do not considerably affect the 2^ nucleon-nucleon 
potential induced by the QCD 9 term, and therefore play a minor role in the calculation of 2^ 
nuclear observables. 


A similar conclusion can be drawn by inserting the value of the gori coupling in the results of Ref. [112] 


where g exchange was considered explicitly. 
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For the phenomenologically interesting coupling gi, all LO relations to baryon mass split¬ 
tings and sigma terms obtain 0(100%) corrections already at next-to-leading order and higher- 
order corrections are even larger. We see no pattern of convergence and conclude that SU{3) 
xPT does not provide a reliable method to extract a value of gi. This coupling plays an im¬ 
portant role in many interesting^observables such as the deuteron EDM and the ^^^Ra Schiff 
moment, such that the lack of a robust relation to the baryon spectrum is unfortunate. In this 
case, the SU{2) xPT extraction discussed in Sect. V C is more reliable but nevertheless suffers 
from a large uncertainty. This uncertainty could be reduced by more precise evaluations of 
the pion mass splitting induced by the quark mass difference. 

As a byproduct of our study, we have obtained expressions for the octet baryon masses 
and mass splittings at N^LO in SU{3) yPT. For diagrams involving octet intermediate state, 
our results reproduce the findings of Ref. m- We also included the effects of the decuplet 
baryons on the mass splittings of the nucleon, S, and S baryons. The N^LO expressions of the 
octet masses and mass splittings depend on several LECs, which cannot be determined purely 
from experimental data. Given the poor convergence/lack of convergence of SU{3) baryon 
xPT |62l |63l I113[ 1114] , it is not clear they can meaningfully be determined from a comparison 
with lattice QCD either. 

In this work we have focused on strong CP violation, but it would be interesting to ex¬ 
tend the study to higher-dimensional GP-violating operators. In many scenarios of beyond- 
Standard-Model (BSM) physics, large nucleon and nuclear EDMs are induced by light-quark 
chromo-electric dipole moments (qCEDMs). However, the sizes of the nucleon EDMs and 
T pion-nucleon couplings ^o,i are poorly known [HI 121] , leading to large uncertainties in the 
analysis of EDM constraints on BSM physics (see for instance Ref. |115j b Just as for the 6 
term, it is possible to derive leading-order relations between qCEDM-induced ^ pion-nucleon 
couplings and baryon mass splittings induced by CP-even quark chromo-magnetic dipole mo¬ 
ments [28|, the chiral partners of the qCEDMs. The baryon mass splittings can be evaluated 
on the lattice providing a method to accurately evaluate go,i- However, the relations have only 
been studied at leading order and they might suffer from large higher-order corrections |116j . 
Finally, a recent evaluation nm of the neutron EDM in SU{3) xPT found a much larger 
dependence on the strange qCEDM than previous studies based on QCD sum rules [H21[nH|. 
As qCEDMs typically scale with the quark mass, this would strongly impact neutron EDM 
constraints on BSM scenarios. However, the analysis of Ref. m is based on leading-order 
SU (3) yPT and, as demonstrated in this work, higher-order corrections might strongly affect 
the results. 
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Appendix A: N^LO corrections to mass splittings of the S and H baryons 


In this appendix we give the expression of the loop contributions to and at 


N"^LO. The counterterm contributions are given in Eq. (72), and the violation to the Coleman- 


Glashow relation in Eq. (74|). The loop functions gi and g 2 are dehned in Eq. 
f 2 {x,y) and fHx.y) in Eqs 


(81) and (82) 


(jTSj), while 

The octet contributions agree with Ref. ca, 


while the decuplet corrections are new results. 


The contribution of the relativistic corrections and of the two-pion couplings bi, 62 and 63 is 
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The decuplet contribution to the S mass splitting is 
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Loop corrections involving the operators bo and bp are 
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The decuplet corrections are 
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Appendix B: N^LO corrections to the octet average masses 


In this Appendix, we give the corrections to the nncleon, S, S and A average masses. In 
order to verify the relations involving the conplings goNiiK and goNAK, in addition to 
A^'^^rriN, A*^^ims and A^^^mA, one needs the octet baryon wave fnnction renormalization, the 
kaon and g meson wave fnnction renormalization, and the corrections to Fk and 

The corrections to the mesons wave fnnction renormalization and decay constants are 


SZk 

SFk 

SF^ 


(1 + T.) < 


64772^2 


32772^2 


Tfl 


16F 

[rris + 2m)L4 + -jy^g + m)L^ 


Fi 
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1287r2Fo2 


(1 + Lk) 


(1 + Lt^) 


F2 
-^0 

3m|- 

647r2F(2 


L^{ms + 2m) + -{2ms + m)L^ 


(1 + Lk) 


3mz 


USti^F^ 


(1 + Lrj) 


8 B 1, 

( (m^ + 2m)L4 + -(m^ + m)L^ 


F2 

-^0 . 

3m|- 
32^^ 


8B 

(1 + Lk) - ( (m* + 2m)Li + ^(2ms + m)L^ ) , 


(Bl) 

(B2) 

(B3) 

(B4) 


where we have introdnce the shorthand Lj = log/r^/mf, for i = {tt, F, r/}, in order to make 
the formnlae in this Appendix more compact. 

Then, we give the corrections to the baryon octet masses and wave fnnction renormalization. 
The resnlt of diagrams [^a) -|^f), which involve octet intermediate states, agree with Ref. 
[65l[75]. The decnplet corrections agree with the resnlts of Ref. [65]. We also agree with [6^ . 
after we expand in A, and set the light qnark mass m to zero, as was done in Ref. |69j . 


• nncleon 


The nncleon wave fnnction renormalization is given in Eq. (85). The loop contribntions to 
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the nucleon mass, including decuplet corrections, are 
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(B 5 ) 


The counterterm operators in give 


(45)^ 


= — ((di — d 2 F ds — F de F dj F d^)ml — { 2 di — 2 dz F d^ — Sde — Ad 7 )mms 
+ (di F d 2 F d^ F 2 d^ + 2 (ig + Adj F 2 dg)ffi^^ . (B 6 ) 


The wave function renormalization of the S held, in the isospin limit, is 


777 "^ 777 ^ 
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The loop contributions to the S mass, including decuplet corrections, are 
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The counterterm operators in give 

2\{ct)j^„ 
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— — ((*^1 T (^2 T T d^ T dg T d7 T d^vfig — (2di — 2dy, — dg — 3d% — Ad'j)TfiTns 
+ (di — d2 + dg — 2d5 + 2d6 + 4d7 + 2d8)m^) . (B9) 


The wave function renormalization of the S field, in the isospin limit, is 
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The loop contributions to the S mass, including decuplet corrections, are 
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The counterterm operators in give 
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A 


The wave function renormalization of the A held, in the isospin limit, is 

(1 + 3L.) + (1 + 3L.) + (1 + 3L,) 

+ 32^2^2 (2/2(mx, A) + 3/2(777,,, A)) . 

(B13) 
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The loop contributions to the A mass, including decuplet corrections, are 
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The counterterm operators in give 
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